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I. MANY-BODY HAMILTONIAN OF GENERIC DIRAC MATERIALS WITH BCUKLING 


In a generic two-dimensional (2D) Dirac material with bcukling, the eigen wave functions are given by Eq. (2) in 
the letter. We divide the wavefunctions into those for two pseudo-layers between which there is a separation d. The 
electrons essentially belong to a double-layer system. Using the density matrix defined in Eq. (3) in the letter, the 
many-body Hamiltonian in the Hartree-Fock approximation (HFA) is given by 


H = ( E o + (°))^ ) Pa,a ( 0 ) 


i ^. () 'y y ^ ^ ^77,77' (c[) i K P G (J ( Q.))Pcr' ,a' (^) 
a,a' q 


E X X V,V' (Pa,a' (“Q ))Pa',a (q) 

a,a' q 


(i) 


where k is the dielectric constant, a, a' = 1 —> ( K , f), 2 —> ( K , J.), 3 —> ( K ', f), 4 —> ( I \', 4_) are the valley-spin indices, 
rj and rf are the valley indices in a and a' respectively, and y x = y 2 = a 2 ,y 3 = = —b 2 . E a is the kinetic energy of 
level a with the Zeeman coupling. The summation with a bar excludes the term of q = 0. The functions H r) , and 
X v r/ , describe the Hartree and Fock interactions between valleys r) and rj ', 
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where J is the Bessel function and 

= a 4 f ntn + b 4 f n _ ln _ 1 +2a 2 b 2 f nn _ 1 e- qd/i , 
Sr,* = (« 4 /n,n + fe 4 /„-i,„-i)e-^ + 2a 2 & 2 /n, n -i, 
with rj being the valley other than rj. The function / is 

fn,m (<?) = e~ q2/2 L n (g 2 /2) L m (q 2 / 2) , 
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with a Laguerre polynomial L n ( L n<0 = 0). 

For a liquid phase at quarter filling of a Landau level the ground state is polarized. The nonzero density matrix 
elements give the distribution of electrons (p ±1 (0)), (p 3 3 (0)) and the coherence (p ± 3 (0)), (p 31 (0)) 


(Pl,l ( 0 )) + <p 3 ,3 ( 0 )) = 1 , 

(Pl,3 (0)) = (P3,l (0)) = \] (Pl,l (0))(P3,3 (0)}> 
where we choose all order parameters to be real. Hence, the energy is given by 
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By using the conditions X K K , = X K , K , X K K = X K , K , and ignoring the constant kinetic energy, we obtain 

e 2 


K,K' ~ ^K' ,Ki ^K,K ~ K',K' 

E = 2~-A(a 2 - b 2 ) 2 d/£ - X KK ( 0 ) + X KK , ( 0)1 <p M (0)) «p M (0)) - l) . 


k£ . 


( 10 ) 


II. MANY-BODY HAMILTONIAN OF SILICENE 


The elements of the density matrix with extra orbital index can be obtained by modifying Eq. (3) in the letter: 


Pa,o;/3,a' (*!) “ Tf E/ e 2(A ’ 1+A ' 2, i Yl c /3,rr',X 2 


( 11 ) 


where a, /? are valley indices, o, o' are orbital indices. Since the orbitals are not conserved the many-body Hamiltonian 
in the HFA is more complex than that in Eq. ©, 

e 2 ~ 
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where E a 0 contains the kinetic energy E a 0 of the orbital o in valley a and a capacitive energy, 
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with the filling factor v. The Hartree interaction is 


h:;Eo 3 ,o 4 (q) = E E E G &,»M 
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where we define a function P m n (q) = e l m n l9 d /^ I 2 anc j coefficients 
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with the angle 6 between vector q and x axis. Function 0 is defined as 

\/min(fc,Z)! (ip) lk ~ 11 L \ k -i\ ( p 2 


The Fock interaction is 
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Once we obtain the HF Hamiltonian, we could obtain the order parameters of the ground state by solving the equation 
of motion of the single Green’s function^. 
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III. PSEUDO-SPIN STIFFNESS OF THE ANISOTROPIC NONLINEAR a MODEL 


We extract the Fock energy functional from the Hamiltonian, 

e 1 2 

H-F = — 'y ' y ] ^oi^Oi,o 3 ,o 2 (*l) (P a , °i;/ 3 ,02 ( — ( l)) Pp,o 3 \a,Oi (q) , 

OL,fi 01,02,03,04 q 


( 20 ) 


where we neglect the kinetic energy, the capacitive energy and the Hartree interaction. This is because the kinetic 
energy is a constant, capacitive energy is very small, and the Hartree term is zero in the slow varying density 
approximation. 

This field theory is only valid when E z = 0. If E z ^ 0, then the ground state is nolonger an easy-plane QHF. The 
ground state will be valley polarized very rapidly, i.e. the E z is very small to polarize the valley pseudo-spin. In the 
finite E z case, the best way is to calculate the microscopic Hamiltonian in the symmetric gauge. 

If we consider the two valleys only in orbital o = N, then the orbital index could be neglected, and 
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1 e 'X'X 
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Then the Lagrangian of the anisotropic Nonlinear a model is obtained by calculating the excitation energy when the 
density matrix is slowly varied^, 


L 


\ Y 


where the normalized pseudo-spin field m is defiend by 

m x (r) + im v (r) = Att£ 2 {p kk , (r)), 

m z (r) = 4t t£ 2 ((p K , K (r)) - {p K , yK , (r))) , 

and the pseudo-spin stiffnesses are given by 
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The excitation energy of a bimeron or an anti-bimeron is 

An r 

SE = -y (Px + P y + Pz) ■ (27) 

Note that this energy which is obtained in the field theory is not identical to a single charged excitation of the electron 
gas. However, the excitation energy of a bimeron-antibimeron pair in the field theory is identical to that in the electron 
gas, which is double of SE, 


6E„ 


87r , 

= Y ( Px 


Pz) ■ 


(28) 


This energy is related to the transport property of the system, which can be measured in a transport experiment. 
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We study the ground states and low-energy excitations of a generic Dirac material with spin-orbit 
coupling and a buckling structure in the presence of a perpendicular magnetic field. The ground 
states can be classified into three types under different conditions: SU(2), easy-plane, and Ising 
quantum Hall ferromagnets. For the SU(2) and the easy-plane quantum Hall ferromagnets there 
are goldstone modes in the collective excitations, while all the modes are gapped in an Ising-type 
ground state. We compare the Ising quantum Hall ferromagnet with that of bilayer graphene and 
present the domain wall solution at finite temperatures. We then specify the phase transitions and 
transport gaps in silicene in Landau levels 0 and 1. The phase diagram strongly depends on the 
magnetic field and the dielectric constant. We note that there exists triple points in the phase 
diagrams in Landau level N = 1 that could be observed in experiments. 


Recently, several graphene-like systems such as sil¬ 
icene, germanene [li] and transition metal dichalcogenides 
(mos 2 ) a have received considerable attention. They 
all have a honeycomb geometry in the xy plane as in 
graphene m , but additionally with a buckled structure 
in the z direction. The buckling is induced by the atoms 
that are heavier than the carbon atoms in graphene. The 
heavy atoms also have more complex electron orbitals. In 
these systems the spin-orbit (SO) coupling is important, 
while the SO coupling is negligible in graphene. Hence, 
the electrons in these systems must be described by a 
massive Dirac equation in which the mass is induced by 
the SO coupling. We generically call these materials as 
buckled Dirac materials. The Brillouin zone is similar 
to that of graphene: a hexagon with two inequivalent 
valleys K and K'. In MoS 2 the T point is important be¬ 
cause the energy near this point is close to that in the 
K , K' valleys in the valence band. For simplicity, in this 
paper, we consider only the K and K' valleys in order to 
compare with graphene. 

In the presence of an applied perpendicular magnetic 
field the electron bands split into a series of Landau levels 
(LLs) near the two valleys. The fractional quantum Hall 
(QH) effect 0 has been studied recently in silicene and 
germanene [f|. The fractal butterflies have also been in¬ 
vestigated theoretically in these systems [j|. Here we re¬ 
port on the ground states and transport properties of the 
symmetry broken states in the integer QH effect regime. 
In bilayer and multilayer graphene )9|, |l0(, earlier theo¬ 
retical works have indicated that the ground states in the 
N ^ 0 LLs are Ising quantum Hall ferromagnets (QHFs), 
since the interlayer Coulomb potential is different from 
the intralayer one. The resulting transport properties 
of bilayer graphene were also observed in an experiment 
Q3- In buckled Dirac materials the buckling divides the 
system into two “pseudo-layers”. Atoms A and B be¬ 
long to different pseudo-layers, respectively. Hence the 
buckling structure makes these monolayer Dirac materi¬ 
als similar to bilayer graphene, i.e., we could observe the 
Ising QHF in these one-atom-layer systems. We discuss 
below the various QH states and the collective modes of 
different QHFs in a few LLs. 


The buckling also couples to an external electric field. 
Without the magnetic field, silicene and germane may 
be converted to topological insulators in an appropriate 
electric field [IH. In the QH regime the phases and trans¬ 
port properties are also much richer and more interesting 
when the electric field is applied. In this work, we will 
specifically discuss how the electric field and the dielec¬ 
tric constant (of different substrates) change the phase 
diagram and control the spin and valley pseudo-spin in 
silicene. These materials are therefore potential candi¬ 
dates for application in spintronics. 

We first consider a generic monolayer Dirac material 
with the SO coupling. The Brillouin zone is in general 
a regular hexagon (as in graphene) with two inequiva¬ 
lent valleys K and K' . The noninteracting Hamiltonian 
around valley y = K , K' in a magnetic field is 

H n = V F {° x P x - WyPy) - ^so^zi (!) 

where 77 = 1 , — 1 for the K and K' valley respectively, v F 
is the Fermi velocity, cr is the Pauli matrix, P = p + eA is 
the canonical momentum and A so describes the SO cou¬ 
pling parameter. The SO strength is also described as the 
mass of the Dirac fermion near each valley. We choose 
the Landau gauge of the vector potential A = (0, Bx, 0). 
The LL energy spectrum is given by E 0 = ±A so and 

E n^o = s S n ( n ) \Aso + 2 {v F h/£) 2 \/ n2 + M, where 
i = sJhfeB is the magnetic length and n is the LL index. 
The eigen wavefunctions in the two valleys are 


if 
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■ V — 


a ^n,X \ 

bh n - i,x' 




K' 

■ \ ~ 


' hh n- 1, X \ 
. a K,X ) 


( 2 ) 


where A' is the guiding center, h n x ( h n<0 = 0) is the LL 
wave function of a two-dimensional electron gas (2DEG) 
in a conventional semiconductor. We define a = \a\ ,b = 
|6|, so the normalization condition is a 2 +b 2 = 1. In MoS 2 
the F point should be included in the Hamiltonian. The 
low-energy effective Hamiltonian is due to the three d- 
orbitals of Mo atoms which are located in the same plane 
fg]- So the MoS 2 is equivalent to a monolayer system 
without buckling. 
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If the geometry of the monolayer Dirac material is ex¬ 
actly the same as graphene then the valley pseudo-spin 
has a SU(2) symmetry in any LL. However, we need to 
consider the buckling when we calculate the Coulomb 
interaction. The ground states of bilayer or multi-layer 
graphene in the N ^ 0 LL are valley pseudo-spin Ising 
QHFs H, HO]. The SU(2) symmetry of valley pseudo¬ 
spin is broken to a Z 2 symmetry because there is a fac¬ 
tor e~ qd difference between the inter-layer and intra-layer 
Coulomb potentials, where q is the momentum and d is 
the distance between two layers. We follow the formal¬ 
ism in SE3 to present a more general classification of 
the QHFs of the ground states in a buckled Dirac mate¬ 
rial. We assume a buckling d in the z direction between 
the two elements of the wavefunction spinors in Eq. ©. 
The buckling divides the wavefunctions into two pseudo¬ 
layers. The density matrix p in the momentum space 
is 


1 


Pa,a' (q) = £ e-^(*i+^ 2 )S Xi 


* ww 2 


X 2 +gJ 2 C cr,X iCt',.Y 2 


(3) 


where a, a' = 1 —> (. K , f), 2 —► (K , /.), 3 —> (A'', t), 4 
(A',|) are the valley-spin indices, the LL degeneracy is 
N,, and the creation and annihilation operators of elec¬ 


trons are c' , c. The average values of the elements of the 
density matrix fully describe the system with the Hamil¬ 
tonian in the Hartree-Fock approximation (HFA) [13]. 

For a quarter filling of a LL which is equivalent to the 
case of three-quarter filling due to the electron-hole sym¬ 
metry, the ground state satisfies {p x 1 (0)) + (p 3 3 (0)) = 1. 
The energy of the liquid phase, without a constant, is 
given by E = 2Q{p lx (0)} ((pip (0)} — l) e 2 / nt, where k 
is dielectric constant, and 


Q (q) — (° 2 — k 2 ) dji — X KK (q) + X K K , (q), (4) 


and Q = Q (q = 0). The first term in Eq. © is the ca¬ 
pacitive energy. The Fock interaction functions X K K (q) 
and X K K , (q) are given by Eq. (3) in [13]. We can clas¬ 
sify the ground states as follows. If Q = 0, the ground 
state is a SU(2) QHF, since (p ±1 (0)} and (p 33 (0)} could 
be of any value to minimize the energy. If Q < 0 , the 
ground state is an Ising QHF. The energy is minimized 
when (p ± 1 (0)) is either 0 or 1. Finally, if Q > 0 the 
ground state is an easy-plane QHF, i.e., the energy is 
minimum only when (p 1 : (0)) = {p 3 3 (0)} = 1/2. 

We define the two-particle Green’s function 


,W,7,7' ( < 1> q'; r l - As) = N p(Pa,a' (~q')> 

-N^(Tp aa , (q,r 1 )p 7) y (-q,r 3 )}, 

where T is the time order operator to study the col¬ 
lective behavior of the system. We employ the gener¬ 
alized random phase approximation (GRPA) to solve 
the equation of motion of the two-particle Green’s 
function. Details can be found in in 0. The 
non-zero lowest-energy collective mode which is given 


by the poles of the retarded Green’s functions is 
c (q 0 ) = 73 \Q [<Pi,i (°)> - (Ps,3 (0))] I , which for 
small q 1 C (q) ~ q 2 . This collective mode, which is sim¬ 
ilar to that of an easy-plane QHF in a double quantum 
well system without tunnelling, is a precess mode of the 
valley pseudo-spin in the xy plane [lj, Hal. In the three 
types of QHFs, when q — > 0 , the collective modes are 
distinguished by their gaps. In a SU(2) QHF, Q = 0, so 
C (q —^ 0 ) — 0. It is a goldstone mode. In an easy-plane 
QHF, (Pip (°)) - (/°3,3 (°)) = 0, so C (q -> 0 ) = 0, which 
means a goldstone mode still exists. In an Ising QHF, 
Q ^ 0 and |(Pi j i (0)} — (p 3 3 (0))| = 1. The goldstone 
mode disappears and all modes are gapped. 

We now consider the two actual materials, silicene and 
germanene, in a perpendicular electric field E z . The elec¬ 
tric field can control the phases and the spin polarization, 
useful for application in spintronics or valley pseudo- 
spintronics. The low-energy noninteracting Hamiltonian, 
in the basis {A t, B f, A B |} is |1T] 

H v = v F ( P x t x - pp y T y ) + r/T z h + dE z r z /2, (5) 

where h = -A so a z - a 0 X R (p y a x - p x a y ) , p = 1 for 
valley K and —1 for the K' valley, r and a are the Pauli 
matrices corresponding to the sublattices and the spin, 
a 0 is the lattice constant, X R is the Rashba SO (RSO) 
coupling and the buckling is d. For silicene (germanene) 
fTsL ITil. these parameters are v F = 5.5 x 10 5 (5.09 x 10 5 ) 
m/s, a 0 = 3.86 (4.06)A, X so = 3.9 (43)meV, X R = 0.7 
(10.7)meV and d = 0.46 (0.66)A. 

The Hamiltonian in Eq. © is more complex than that 
in Eq. ( © since the electric field and the RSO coupling 
are included. In the QH the wave function in valley a and 

orbital o is (c“ 1 h 0+ai c“ 2 h 0+a2 c“ 3 /i 0+a;j c“ 4 h 0+Q J , 

with the normalization condition X)i=i | c o iT = A where 
K 1 = I< 4 = 0, K, 2 = -1,K 3 = 1; K[ = -1,A' = A' = 
0, K 4 = 1. Because of the RSO coupling the eigenvec¬ 
tors are not spin polarized. We introduce another de¬ 
gree of freedom, the orbital, to replace the spin. We find 
that without Zeeman coupling the energies of the orbitals 
o = N, N — 1 are close to each other in the LL N. The 
concept of the orbital degree of freedom here is similar 
to that in the N = 0 LL in bilayer graphene. The RSO 
interaction couples different spins in a valley. For zero 
RSO coupling the orbital degree of freedom is identical to 
spin. In reality, |c^-_ 1 ;1 | , |c^_ 12 | , |c^ 3 | , |c^ 4 | < 10“ 4 , 
and so approximately, the orbital N is associated with 
spin up and orbital N — 1 is associated with spin down. 
Note that the coefficients c“ 4 not only depend on the 
magnetic field but also on the electric field. 

We neglect the LL mixing since the LL gap is large 
(E n =o ~ 0 ,E n _ ±1 ~ ±60 meV for silicene). The 
density matrix and the many-body Hamiltonian in the 
HFA are given in 0- We define the Green’s function 

G a ,o-,p,o> ( x ,X'x) = -{'TCa, 0 ,x( T ) 4, 0 ', X ’ (0))’ with a 

relation to the density matrix at zero temperature, 
G a ,o-p,o' (q± = °“) = (Pp,o'-, a ,o (q)>; Solving the equa¬ 
tion of motion of the Green’s function we could obtain 
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FIG. 1: (Color online) (a) The occupation of the states |1) 
and 13), and the valley polarization ( p z ) at v = — 1. (b) Order 
parameters around the phase transition region at v = 0. 

the ground states of the system 0]. In what follows, we 
define the valley-orbital (or say, valley-spin if X R = 0) 
indices as ( K,N ) —>• 1, (K, N — 1) —► 2 ,(K',N) —> 
3, (KN — 1) —> 4, for simplicity. 

Due to the electron-hole symmetry, we consider only 
the quarter- and half-filled LLs. The system can be 
described in the (pseudo-)spin language. The valley 
pseudo-spin field in orbital o is defined by p o x + W 0 . v = 

(Pk,o-K',o) iPo,z = (PK,o\K,o) ~ (pK',o;K',o) > an d P = 
S 0 Po- We could approximately associate the orbital 
with the real spin: o = N associated with spin up and 
o = N — 1 associated with spin down. Hence, we de¬ 
fine the spin field: S a x + iS ay = (p„,jv ; a,iv-i) , S a,z = 

( Pa,N-,a,N ) — (Pa,N-l-,a,N-l) i an ^ S = Y) a 

Filling factor v = —1 : When E z = 0, the ground state 
is an easy-plane QHF in valley pseudo-spin. It can also 
be obtained by the classification parameter Q in Eq. 

Q > 0 if we approximate X R = 0. The ground state lies 
on the easy-plane QHF regime. It is the symmetric state 
of the two valleys IGS 1 ) = (|1) + |3)) /y/2. When the 
electric field increases, there is a bias between the two 
states |1) and |3) so that the ground state is a bonding 
state \GS) = a 1 11) + a 3 |3). When the electric field is 
strong E, > 0.08 mV/nm, the bias is large enough to 
polarize the valley pseudo-spin. The order parameters in 
the phase transition are shown in Fig. |T| (a). 

For an easy-plane QHF, the charged excitation is a 
bimeron or an anti-bimeron described in the anisotropic 
nonlinear a model [53, [22j]. At E z = 0 the La- 
grangian of this model can be written as L = 
5 y z Pn > w h ere we define the normalized 

field m =47r£ 2 p. The pseudo-spin stiffnesses p M are de¬ 
fined in Eqs. (25) and (26) in 0- The excitation en¬ 
ergy of a bimeron-antibimeron pair is then given by 
SE pa i r = 8n (p x + p y + Pz) / 3 [Hi 01 • For B = 10T, 
E z = 0, k = 1, the excitation energy of a bimeron- 
antibimeron pair is 39.2 meV. In comparison, the excita¬ 
tion energy of an electron-hole pair is 156 meV. Hence, 
the transport gap is due to the bimeron-antibimeron pair. 

Filling factor v = 0 : For E z = 0, |1) and |3) are 
fully occupied, the ground state is spin polarized and 
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FIG. 2: (Color online) The excitation energy contour of the 
quasi-particle around the triple point at v = 3, k = 5. The 
dashed lines are located at the phase transitions. The red dot 
is the triple point. 


is valley unpolarized. The ground state is stable when 
E z < 219.05 mV/nm. In the region E z £ [219.05, 219.25] 
mV/nm, the coherence (p x 4 ) between |1) and |4) arises. 
In Fig. (T] (b), we show the phase transition region where 
(p 1 1 ) and (P 44 ) are gradually changed. The spin of the 
system is also controlled unpolarized gradually by the 
electric field. When E z > 219.25 mV/nm the system 
is spin unpolarized but valley pseudo-spin polarized, all 
electrons are in valley K'. 

In the N > 0 LLs the nature of the broken symmetry 
states are different from those in the N = 0 LL. For 
N > 0 the ground state of a LL is an Ising QHF, which 
is similar to bilayer graphene 00 - For simplicity and 
without loss of generality, we only study the filling factors 
v = 3,4, since the LL mixing is important in higher LLs. 

Filling factor v = 3 : For E z = 0, we assume that 
X R = 0 so that Q < 0. Hence the ground state is an Ising 
QHF with a valley Z 2 symmetry, which is also supported 
by our numerical calculation including X R . The SU(2) 
valley symmetry is broken to a Z 2 symmetry. Figure [ 2 ] 
shows the phase diagram in a finite electric field for k = 5. 
For X R = 0 the SU(2) spin symmetry is also broken to 
a Z 2 symmetry at the phase transition between Phase 
I (or Phase III) and Phase II. These symmetry broken 
states are all induced by the small buckling geometry. 
Note that the valley and spin can also be controlled by 
the electric or the magnetic field. At the two sides of 
the phase transition line in Fig. [2] the valley or spin is 
reversed. 

Interestingly, for B = 15.5 T and E z = 1180 mV/mn, 
there is a triple point (the red dot in Fig. [2| in the phase 
diagram. This triple point occurs only when n > 3. 
When the dielectric constant is very large the electron gas 
is close to a noninteracting system and the triple point 
would disappear. At a finite electric field the phase di¬ 
agram is also changed by the dielectric constant n. For 
k = 1 the phase III disappears and only other two phases 
survive when E z < 2000 mV/nm. 

The Ising QHF here is similar to that in bilayer or mul¬ 
tilayer graphene 0 , but is different from the Ising QHF 
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FIG. 3: (Color online) The phase diagrams at v = 4 and (a) 
k — 1, and (b) k = 4. 


system with different LLs [22$ in semiconductor quan¬ 
tum wells. Hence, the lowest charged excitation may be 
a skyrmion around E z = 0. However, the skyrmion in 
this system must be calculated numerically with a mi¬ 
croscopic Hamiltonian in the symmetric gauge. 

There is no domain wall at zero temperature. At fi¬ 
nite temperature the domain wall could exist to lower 
the free energy of the system [ 2 J], when we consider 
the wall entropy. Below a critical temperature Tc, do¬ 
main walls provide ID channels carrying extra charges 
(electron-hole pairs) to dissipate the transport charge of 
the 2DEG, when the domain walls are dense enough to 
overlap. So the resistance spike in R xx appears. Above 
Tc, the domain wall will be infinitely long and expand to 
the sample perimeter. The charge in the domain wall can 
not dissipate the transport electrons any more and hence 
the resistance spike disappears. Following the study of 
the domain wall in a graphene bilayer 0 we obtain the 
kink domain wall of the valley pseudo-spin at E z = 0: 
m x (r) = sin 9 (r), m y = 0 and m~ (r) = cos 9 (r) with 


9 = 2 arctan exp 


^2(K z ~K ± )/p s x 


( 6 ) 


where we approximate p s = p x ~ p z , and define Kp = 

■^■n’,n,n,n (0)/(^ 7r ) an d k-z = X n ’n,n,n (0) /(^ 7r )- The 
excitation energy per unit length of the domain wall is 
then given by SE = 2yQ ( K z — K ± ) p s . 

We also present the quasi-particle (QP) excitation en¬ 
ergy Eq P around the triple point in Fig. [2] Experimen¬ 
tally, the phase transitions between different spins (Phase 
I or III to Phase II) can not only be observed in a NMR 
experiment, but also in a transport measurement. The 
phase transition between different valleys (Phase I to III) 
may be observed in a transport or a circular light absorp¬ 
tion experiment [25[. At the phase transitions a resistiv¬ 
ity spike may be observed due to the existence of the 
domain wall at finite temperature, which has been re¬ 
ported in a LL Isin g QH F [23j, [24j and in the Ising QHF 
of bilayer graphene [HJ, 111 ]. 

Filling factor v = 4 : For k = 1, the phase diagram in 
a magnetic field is shown in Fig. [3] (a). The triple point 
appears when k > 3. We show an example of the triple 
point which is marked as a red dot located at B = 9.2T, 


FIG. 4: (Color online) The collective modes for (a) the easy- 
plane QHF ground state at v = —1 and (b) the Ising QHF 
ground state at v = 3, when B = 10T, k = 1 and E z = 0. 
The small q region is given as inset. 


E z = 635mV/nm for k = 4 in Fig. [3] (b). Phase C sets 
in when B < 9.2T, since the kinetic energy contributes 
more and the Coulomb interaction e 2 /«£ is decreased by 
the large dielectric constant. 

As we discussed above, the easy-plane QHFs have a 
goldstone mode while all the modes of Ising QHFs are 
gapped. In silicene we only show the collective modes at 
v = —1 and 3 in Fig. [J] for simplicity. Note that, in the 
region E, £ [219.05,219.25] mV/nrn, the ground state is 
a bonding state with a goldstone mode at v = 0. 

To summarize, we classify the ground states of a 
generic buckled Dirac material in a magnetic field into 
three types of QHFs. The low-energy collective modes of 
the three QHFs are given analytically in the GRPA. A 
goldstone mode exists in the SU(2) and in the easy-plane 
QHFs, but not in the Ising QHF. We then focus on a 
real material, viz., silicene. Without an electric field we 
note that in silicene the magnetic field is able to change 
the coefficients cf i in the wavefunctions. However, in a 
very small magnetic field (B -C 0.01T) the ground state 
becomes an easy-plane QHF at v = 3. In such a low mag¬ 
netic field the QH effect can not be realized and the LL 
mixing is not negligible. If the SO coupling can be tuned 
then the coefficients of wave functions can be modified. 
For B = 10T, v = 3, k = 1, we find that when X so > 750 
meV, the ground state at v = 3 is an easy-plane QHF. 
For about X so = 750 meV, the ground state is a SU(2) 
QHF. This is also true for germanene. If we could effi¬ 
ciently tune the wave function we may realize the phase 
transition between different QHFs. Experimentally, this 
transition is observable: the domain wall induced resis¬ 
tivity spike occurs only in an Ising QHF. The phase di¬ 
agrams and transport properties in the N = 0,1 LLs 
in silicene also depend on the magnetic field and the di¬ 
electric constant which dramatically change the Coulomb 
interaction. We have shown the triple points in Figs. [ 2 ] 
and[3| The SU(2) symmetry of the spin and valley are 
broken by the electric held and the buckling structure. 
The phase transitions may indeed be observed in NMR, 
transport or light absorption experiments. 
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